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Abstract 

The shelling of crystals is concerned with counting the number of atoms on spherical shells of a given radius and a 
fixed centre. Its straight-forward generalization to quasicrystals, the so-called central shelling, leads to non-universal 
answers. As one way to cope with this situation, we consider shelling averages over all quasicrystal points. We express 
the averaged shelling numbers in terms of the autocorrelation coefficients and give explicit results for the usual suspects, 
both perfect and random. 
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1. Introduction 

One characteristic feature of a crystal is the num- 
ber of atoms on shells of radius r around an arbitrary, 
but fixed centre, e.g. around one fixed atom. In the 
simplest case, one thus considers a lattice, such as 
the square lattice Z 2 in the plane, and determines the 
number of lattice points on circles of radius r. In the 
square lattice case, only squared radii r 2 = to 2 + n 2 
with to, n € Z are possible for obvious geometric rea- 
sons, compare Fig. |l|. 

For this example, the answer is well known from 
number theory [|| because it essentially boils down 
to counting the number of ways that r 2 can be writ- 
ten as the sum of two squares. The result reads as 
follows. Let r 2 = M be an integer, and a(r) the 
number of points of Z 2 on a circle of radius r around 
the origin. Then, a(r) = Aa(M) where a(M) is a 
multiplicative arithmetic function, i.e. a(l) = 1 and 
a{MN) = a{M)a{N) for coprime M, N. So, it suf- 
fices to know a{M) for M a prime power, i.e. M = jr. 
The explicit result is 




if p = 2 

if p = 1 mod 4 

if p = 3 mod 4 and I odd 

if p = 3 mod 4 and I even 



(1) 




Fig. 1. Shelling of the square lattice. The circles centred 
at the lattice point denoted by the thick dot have radii 
r 2 = 1, 2, 4, 5, 8, 9, 10, 13, 16, the corresponding shelling 
numbers are a(f) = 4, 4, 4, 8, 4, 4, 8, 8, 4, compare Eq. (|l]). 



see I! for further examples and a systematic exposi- 
tion by means of lattice theta functions. 

Now, if one tries to extend this approach to qua- 
sicrystals, one immediate problem arises: due to the 
lack of periodicity, there are no "natural" centres any 
more, except, perhaps, the centre in a pattern with 
exact maximal symmetry — if this exists! There are 
several attempts to find analogues of Eq. (p| ) , and ex- 
plicit results can be found in jliL |l2, 13 l(J|j|. How- 
ever, an alternative approach is desirable in which 



There are many explicit formulas known for lattices, one does not ask for the shelling numbers around a 
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fixed centre, but for the averaged shelling numbers 
around any possible centre e.g. around any cen- 
tre that is itself in the point set. This is what we 
want to continue to analyze here, in analogy to earlier 
work on the coordination numbers of crystals versus 
quasicrystals [|. Note that the ordinary and the 
averaged shelling problems coincide for lattices. 

2. A formula for averaged shelling 

It is intuitively clear how to define the averaged 
shelling number for a point set A: if we fix a radius r 
and a centred ball B s (0) of radius s > r, we inspect 
any point inside the large ball, determine the shelling 
number around it, and average over these possibili- 
ties. This is a radially averaged autocorrelation coef- 
ficient. The problem with it is that there might not 
be a limit as s — > oo, and this is indeed a serious 
problem in general. However, cut-and-project sets, 
or model sets as we want to call them, have the nice 
property that these averages are guaranteed to ex- 
ist. What is more, each finite patch or configuration 
within a regular, generic model set (see |H| for the 
precise conditions) occurs with a uniform frequency. 
Thus, we can evaluate the averaged shelling number 
by means of the existing autocorrelation coefficients. 
Let us describe how this works. 

Let A = A(f2) be a regular model set [|ll] with 
window il. The possible distances between points 
are all of the form r — \y\ with y £ A — A where 
— x 2 | x 1 ,x 2 G A}. We find the formula 



cr(r) 



E Kf) 
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Here, for y e A — A, v(y) is the autocorrelation coef- 
ficient, normalized per point rather than per volume: 
1 
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where A s = A fl B s (0), \A S \ is the number of points 
in A s , (.)* is the *-map (or li ft) from physical to 
internal space as described in uM, and Xq denotes 
the characteristic function of the window. Note that 
the last step of (^) is correct because Weyl's Theorem 
applies to regular model sets [|l5| . 

3. Result in one dimension 

Let us explain the typical situation in one dimen- 
sion with a representative example. Consider the 
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Fig. 2. Averaged shelling function a(r) for the silver 
mean model set, plotted against the distance r G Z[\/2] in 
physical space (left) and against its counterpart r int = |r'| 
in internal space (right). 



silver mean substitution rule a i— > aba, b i— > a. This 
gives a semi-infinite fixed point of the form 

\abaaabaabaabaaabaabaaabaabaaab ... (4) 

If a codes an interval of length 1 + \/2 and b one of 
length 1, this gives a point set, starting at 0, which 
is the positive part of the silver mean model set 



A 



G Z[V2] 



x'e 


■ \/2 \f2- 




. 2~' ~~2 . 



} (5) 



where (.)' is the algebraic conjugation, defined by 
\/2 i ^ — v2, and represents the *-map in this case. 

A— A is a subset of Z[\/2], and itself another model 
set, this time with window [— y/2, y/2]. So, applying 
Eq. (||), the autocorrelation coefficient for y € A— A 
can be written as v{y) — f(y') with 



m = 



if|y'l<\/2 







otherwise. 



(6) 



The possible shelling radii are then the non-negative 
numbers in A — A. So, the shelling function is simply 
ct(0) = 1 and, since v(y) = v{-y), a(y) = 2f(y') for 
any positive y E A — A. Now, the averaged shelling 
can be calculated, and the smallest positive radii are 
1, 1 + V2, 2 + a/2, 2 + 2\/2, 3 + 2V2 and so on. The 
result is plotted in Fig. |[ It is clearly seen how the 
inherent structure of the averaged shelling is hidden 
in physical space but clearly visible in internal space. 

Let us briefly mention that one can also calculate 
the averaged shelling for the corresponding random 
tiling, as the autocorrelation is explicitly known in 
one dimension, see || for details. 

4. Rhombic Penrose tiling 

In the plane, the most studied non-periodic tiling 
is the classic Penrose tiling in its rhombic version, 
see the left part of Fig. |j| We consider a realization 
where the edges have unit length. The set of ver- 
tex points constitutes a model set, however one with 
four different translation classes and hence also four 
different windows B . Although it is still possible to 



2 




Fig. 3. Patches of the rhombic Penrose (left) and the 
Ammann-Beenker tiling (right). For the latter, two sim- 
pleton flips in the shaded hexagons are indicated. The 
thin lines correspond to the patch of the perfect tiling. 

use Eq. (^), it is computationally a little easier to de- 
termine all finite patches of a certain radius together 
with their frequencies and to calculate the exact aver- 
aged shelling from that. This is done by a refinement 
of the window method as described in Q . The result 
is summarized in Table [jj. 

Note that the squares of possible radii, r 2 , are 
totally positive numbers (i.e. both r 2 and its alge- 
braic conjugate are positive) in the ring Z[r] where 
r = (1 + v / 5)/2 is the golden number. This can be 
understood from the underlying shelling problem of 
the cyclotomic ring Z[£] with £ = exp(27ri/5) and 
its number-theoretic ramifications, see || for details 
and a closed formula for the central shelling of Z[£]. 

5. Octagonal tiling: perfect versus random 

Finally, let us consider the vertex set of the 
well-known Ammann-Beenker or octagonal tiling as 
shown in the right half of Fig. ||. As above, edges 
are fixed to have unit length. The vertex points 
form again a simple model set, obtained (up to scale) 
from the four-dimensional lattice Z 4 with a regular 
octagon as window in the 2D internal space S. To 
obtain a good numerical approximation for a larger 
number of radii in this case, we have used a large pe- 
riodic approximant (with 47321 vertices in the unit 
cell) and determined the averages explicitly. For 
small radii, the agreement with the exact formula 
based upon the window technique is 6 to 8 digits. 
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Fig. 4. Averaged shelling function cr p (r) for the perfect 
Ammann-Beenker tiling, plotted against the radius r in 
physical space (left) and against its counterpart 
internal space (right). 



Table 1. Exact results for the first 50 averaged shelling 
numbers er(r) of the rhombic Penrose tiling. 
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On the left of Fig. |[ the averaged shelling is 
shown against the radius in physical space. The 
result looks as erratic as in the ID case, except 
for the square-root growth of the averaged shelling 
numbers at small distances. On the right, plot- 
ting the shelling against the internal radius shows 
the inherent structure. The points lie on curves 
that start at multiples of eight on the ordinate and 
glide down convexly towards the abscissa which is 
met at a point between = 1 + V2 ~ 2.4142 and 
r 2 = (4 + 2v / 2) 1 / 2 ~ 2.6131. These radii correspond 
to the inradius and circumradius of the window of 
the set A — A, which is an octagon of sidelength two. 

Alternatively, one may consider the corresponding 
random tiling. It can be obtained, starting from the 
perfect pattern, by simpleton thermalization, i.e., by 
repeating random simpleton flips as shown in Fig. |3| 
for sufficiently many times (some 1000 flips per ver- 
tex of the starting patch as a rule of thumb) . Doing 
this for a periodic approximant with 8119 vertices 
in the unit cell and numerically evaluating the aver- 
aged shelling numbers for the resulting random tiling 
gives the other entries of Table 0. Note that a real 
ensemble average, which would give the correct (and 
exact) average, is practically impossible, so we have 
to rely on the self-averaging nature of the model. 
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Table 2. Averaged shelling function a p (r) for perfect 
and <J s {r) for stochastic Ammann-Beenker tilings. The 
possible shell radii of the perfect tiling are a proper subset 
of those of the random tiling. 
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6. Open problems 

Illustrated by explicit results for some of the most 
common examples, we have demonstrated what the 
averaged shelling for quasicrystals means and how 
the corresponding results look like. However, closed 
expressions seem difficult already for planar exam- 
ples, and their real implications are still to be un- 
derstood. It is clear that several interesting number- 
theoretic connections are lurking in the back, at least 
if there is any reasonable generalization of the notion 
of theta functions. Some further hints at this can be 
found in @], but additional work is needed before a 
more definitive statement can be made. 

From the computational point of view, one can use 
the above approach and complete the exercise for the 
other standard examples, in particular for the icosa- 
hedral cases in three-dimensional space and for the 
Elser-Sloane quasicrystal M. One interesting aspect 



certainly is the possibility to distinguish perfect from 
random order this way, as is apparent from the octag- 
onal tilings treated above. This is to be compared 
with the Fourier space approach || [l(J and could 
give rise to improved statements about the validity 
of perfect versus random tiling models. 
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